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We study the problem of controlling a general complex network towards an assigned synchronous 
evolution, by means of a pinning control strategy. We define the pinning-controllability of the 
network in terms of the spectral properties of an extended network topology. The roles of the 
control and coupling gains as well as of the number of pinned nodes are also discussed. 



I. INTRODUCTION 

The need of regulating the behavior of large ensembles 
of interacting units is a common feature of many physical, 
social and biological networks. For instance, many reg- 
ulatory mechanisms have been uncovered in the context 
of biological, physiological and cellular processes, which 
are fundamental to guarantee the correct functioning of 
the whole network dj. Examples include, in physiol- 
ogy, the synchronous beat of heart cells whose rhythm 
is generated by pacemaker cells situated at the sinoatrial 
node , the control of the respiratory rhythm played by 
synaptically coupled pacemaker neurons in the medulla 
[3| and the generation of rhythmic pacemaker currents 
by networks of interstitial cells of Cajal in the gastroin- 
testinal tract of mammals In social networks, opinion 
dynamics are often driven by key-individuals termed as 
opinion leaders Q. Other relevant examples, in biology, 
are the mechanisms through which cell cycles are con- 
trolled [|| and synchronized 

Understanding the fundamental nature of such regu- 
latory mechanisms is therefore of utmost importance in 
Physics and Applied Science. Pinning control has been 
proposed in the literature [1, Q as a fitting model to 
provide an insight into the regulatory mechanisms to 
control lattices and networks of coupled dynamical sys- 
tems. The general idea behind pinning control is a self- 
feedback action (over a given reference evolution) , acting 
on a limited subset of the dynamical systems placed at 
the network nodes. These nodes, also termed as refer- 
ence sites or pinned sites, play the role of network lead- 
ers/pacemakers. Specifically, a direct control action is 
active only on these nodes and is propagated to the rest 
of the network through the coupling among the vertices. 

In this paper, we define the new concept of pinning- 
controllability (different from synchronizability) and de- 
rive a new quantity to assess this property for physi- 
cal networks of interest. These findings have immediate 
theoretical and experimental relevance. By evaluating 
the pinning-controllability of a given network, one could 
decide the effectiveness of a pinning control strategy in 
terms of the strength of the required control action, the 
number of nodes to be pinned and the effects that a given 
topology of feedback connections can have on the entire 
network. 

The rest of the paper is outlined as follows. In Sec. 



II we introduce the pinning control formalism, apt to 
describing dynamical complex networks subject to a de- 
centralized control action. In Sec. Ill a new methodology 
(based on the Master Stability Function approach 
is introduced to define the pinning-controllability of a 
given complex network of interest. Finally, a validation 
is presented in Sec. IV, where pinning control of scale- 
free complex networks of chaotic oscillators is studied 
through numerical simulations. 



II. PINNING CONTROL: AN OVERVIEW 

In recent years, synchronization of complex networks 
of coupled oscillators has been the subject of intensive 
research activity within the scientific community. A com- 
mon assumption in such literature is that all the dynam- 
ical systems at the network nodes are identical, while the 
problem of synchronization of networks of non-identical 
oscillators has received much less attention. In particu- 
lar, the Master Stability Function (MSF) approach, in- 
troduced in has been successfully applied to syn- 
chronization of complex networks in a wide variety of 
situations (see e.g. QiJ HH, EH EH, to name a few con- 
tributions). Note that the MSF approach is valid under 
the hypothesis that all the vector fields at the network 
dynamical nodes are identical. 

In this paper, we study networks in which two different 
layers of dynamical nodes coexist: the uncontrolled sites 
and the reference (controlled) ones. In particular, the 
latter plays the role of leading the whole network towards 
a given (desired) reference evolution. 

Generally, we assume the controlled complex network 
to be described by the following set of equations: 



dxi 
~~dt 



N 



f(xi) + o-y^^Aij{h(xj)-h(xi)) + ay^^6(i-Ck)ui, 

3 = 1 k=l 

i = 1, N, representing the behavior of N identical dy- 
namical systems coupled through the network edges. 

The first term on the right side of ([TJ) describes the 
state dynamics of the oscillator at each node, {xi(t), i — 
1,...,N}, via the nonlinear vector field f(xi); the sec- 
ond term represents the coupling among pairs of con- 
nected oscillators, through a generic output function 
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h{xi), where the coupling gain a represents the over- 
all strength of the interaction. Information about the 
weighed network topology is contained in the coupling 
matrix A, whose entries Aij, are zero if node i is not 
connected to node j ^ i, but are positive if there is a 
direct influence from node i to node j, with | Aij | giv- 
ing a measure of the strength of the interaction. In what 
follows, we assume the matrix A to be irreducible (i.e. 
the associated digraph is globally connected). The last 
term on the right-hand side of ([1]) represents the pinning 
control action. This term is present only for n = pN 
(usually pinned nodes in the network, identified 

by the set C = {c\, C2, Cn}. As commonly assumed in 
pinning control schemes, such nodes play the role of lead- 
ing the others towards some desired reference evolution, 
say s(t). Note that the control input Ui = ip(xi,s,t) has 
a direct influence only on the nodes belonging to the set 
of the reference sites C . 

In what follows, we consider two different strategies for 
choosing the pinned nodes: (i) Random pinning: The n 
pinned nodes are randomly selected with uniform prob- 
ability from the set of all the network vertices, (ii) Se- 
lective pinning: The n pinned nodes are firstly sorted 
according to a certain property at the network vertices 
(for instance, the vertex degree or betweenness central- 
ity), then the nodes to be pinned are chosen in that par- 
ticular order. 

Hereafter, following 0, we choose the control input 
Ui to be generated by a simple state-feedback law with 
respect to the reference evolution s(t), which is assumed 
to satisfy 4| = f(s). Thus, we set Ui — Ki(s — Xi) at 
every pinned node Cj, i = 1, n, where Hi is the control 
gain acting on node Cj. 

Eq. fTJ) can be rewritten as 



dxi 
~dt 



N 



f(xi) - a}^Cijh(xj) -a)^6(i- c^n^Xi - s), 
j=i k=i 

(2) 

i = 1, N, where the elements of the Laplacian matrix 
L, are as follows: Cij = — Aij if i ^ j and La ~ J2j^i ^ij 
Vi. Denote by {Xi = A£ + the set of eigenvalues 
of C and assume they are ordered in such a way that 

xi<\ r 2 <---<\ r N . 

Note that in the specific case where C = 0, i.e. no 
reference sites are present in the network, it is possible 
to define the synchronization manifold, S := {xi = x-i = 
... — xn}, which is an invariant set for system ^ jlQj |. In 
this case, according to the choice of the coupling gain a, 
all the oscillators in the network may settle over a solu- 
tion belonging to S (not known a priori) and the network 
becomes synchronized. Interestingly, in [lfj, lllj , the sta- 
bility of the manifold S was related to R N = and 

M N = maxj Xj , suggesting an evaluation of the synchro- 
nizability of a generic network topology in terms of its 
spectrum. In particular, a network is said to be more (or 
less) synchronizable according to the width of the range 
of values of the coupling gain a for which synchronization 



is attained. 

A different problem consists, instead, in driving a net- 
work of coupled dynamical systems towards a desired 
evolution determined a priori, by means of an appro- 
priate control action. This is the case where C ^ 0, indi- 
cating the presence of some closed-loop equations within 
the set C in (J2]). Note that when C ^ 0, the evolution 
xi(t) = x 2 (t) = ... = xpf(t) = s(t) becomes the only 
admissible solution in S for the controlled network de- 
scribed by ([2]). Thus, it is important to study the stabil- 
ity of such a solution, which corresponds in this context 
to the model reference chosen for the network. 

It is worth noting here that the Master Stability Func- 
tion approach cannot be applied directly to network 
because of the presence of inhomogeneous dynamics at 
the controlled/uncontrolled nodes. In the next section 
we will provide an extension of the MSF approach to the 
case of dynamical networks subject to pinning control, 
described by the set of equations ((2]). 



III. EVALUATING THE 
PINNING-CONTROLLABILITY OF COMPLEX 
NETWORKS 

In this section, we propose for the first time that the 
approach presented in [I0( can be extended to define and 
assess the pinning- controllability of a given network of in- 
terest. In particular, given a network described by ([2]), 
we define the network pinning-controllability in terms of 
the values of the coupling gain a and the control gains 
Ki, i = 1, n, in ([2]), needed in order for the network to 
achieve the desired evolution 4f = /(s). We shall seek to 
assess the pinning-controllability of network @ by con- 
sidering an extended network of iV+ 1 dynamical systems 
yi, where yi = Xi for i — 1,2, .., N and yN+i = s. In so 
doing, we assume that the desired common evolution s 
is given by the state evolution of an extra virtual vertex, 
which is added to the original network. Then, we can 
rewrite ([2]) as 



dyi 
dt 



N+l 



i = l,2,...,JV+l, (3) 



where Ai — {M-ij} is an [N + l)-dimensional square 
matrix, with Mij = 
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and & = YJl=i - c fe)- 

Note that M. is an asymmetric zero row-sum matrix, 
with positive values along the main diagonal. At first, let 
us assume the matrix A to be symmetric, which ensures 
the matrix M. to be diagonalizable (23|. Then, let {/i; = 
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Mi + 3 'Mil be the eigenvalues of A4 and assume they are 
ordered in such a way that /LtJ < /z£ < • • • < /i^- By 
graph theory and linear algebra, we have [i\ > Vi [12l ]. 
Moreover, /ii is the only null eigenvalue of Ai. 

Now, the arguments of the Master Stability Function 
approach [l(| immediately apply to (O . In particular, by 
making use of the MSF, we can investigate the stability 
of the reference evolution, y\{t) = 2/2 (i) = ••• = UN{t) = 
VN+i(t), which can be evaluated in terms of the stability 
of N + 1 independent blocks in the parameter a = afii , 
i = 1, N + 1. Once the functional forms of / and h 
have been assigned, the MSF associates to each value of 
the complex parameter a = a r +ja l , the maximum Lya- 
punov exponent of the system in each block. It can be 
shown that for a large class of systems (in terms of the 
dynamic function / and the output function h), there 
exists a bounded zone of the complex plane centered on 
the real axis, for which the MSF is negative. Thus, the 
condition to be satisfied, in order to guarantee the sta- 
bility of the desired common solution s(t), is that all the 
afii, i = 2, N + 1, belong to the bounded reg ion of the 
complex plane where the MSF is negative [lCj . Namely, 
for a given form of the matrix Ai, there is typically a 
finite range of values of a, say E, for which the stability 
of the synchronous state can be achieved. 

Hence, the spectrum of Ai affects the stability of the 
synchronization manifold. In this paper, we focus on 
the case that Ai has a real spectrum. Then, the lower 
the eigenratio R N+1 = ^|±I, the larger the E [l(J. 
More generally, when the Laplacian matrix has a com- 
plex spectrum, it has been proposed that the width of 
E can be evaluated in terms of both R N+1 — + 1 and 

M N+1 = maxj fij 11. Namely, the condition on the 
maximum imaginary part of the spectrum gives informa- 
tion about the spread of the eigenvalues along the direc- 
tion of the imaginary axis; in the limit of M N+1 — ► 0, 
the entire spectrum tends to the real axis and the (best) 
condition is recovered that E depends only on R N+1 . 
Furthermore, as explained in [13] , the same approach 
can be extended also to the case of Ai being a non- 
diagonalizable matrix, when the condition is satisfied 
that the network embeds at least an oriented spanning 
tree (H) . 

Now, in the case often considered in the literature 
of all control gains being the same, i.e. k = «i = 
k% = ... = re,„, by rewriting the system in ([3]) as 
%t = f(Vi) - M^hfa), for 1 = 1,...,N+ 1, 

where Ai' = ^Ai, and by noticing that the two matrices 
Ai and Ai' are characterized by the same eigenratios, 
we obtain, for every choice of ct, a similar condition on 
the interval of values of re, say /C, for which the reference 
evolution is stable (this is true as long as the maximum 
imaginary part M N+l remains negligible). 

Note that, by taking this approach, we have succeeded 
in decoupling the dynamical properties of the open-loop 
network (in terms of / and h) from the factors encoded 
in the matrix Ai. Specifically: (i) the structural prop- 
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FIG. 1: Left pictures: Combinatorial Laplacian (symmetric 
topology), (a) Behavior of the eigenratio of the 

Laplacian spectrum as function of the control gain re for a 
BA network, with average degree (k) = 2, p = 0.1, a = 0.30. 

(b) Control error at regime E as function of the control gain 
re under the same conditions as in the upper plot. 

Right pictures: Normalized Laplacian (asymmetric topology). 

(c) R N+1 vs. re for a BA network, with average degree {k} — 2, 
p = 0.1, a = 2.8. (d) E vs. re under the same conditions as in 
the upper plot. In (a) and (c) the horizontal continuous lines 
represent the eigenratio R of the corresponding uncontrolled 
networks. 



erties of the network, in terms of the original network 
topology and the choice of the set of the controlled nodes 
C, and (ii) the choice for each node in C, ci, C2, c n of 
the associated control gain k\, K2, K n . 

Thus, in analogy to the concept of network synchroniz- 
ability defined in the literature [l(| [3 EH ; we propose 
to give a definition of the network pinning- controllability 
in terms of the widths of the ranges E and /C, for which 
the reference evolution yi(t) = y%(t) = ... = yN{t) = s(t) 
is stable. Specifically, the lower the R N+1 and M N+1 
are, the more the network is pinning-controllable (note 
that this definition is independent of the choice of the 
functions / and h). Also, according to this definition, 
it becomes possible to act on both the network topology 
as well as the choice of the nodes to pin and their con- 
trol gains, in order to vary (and eventually improve) the 
network pinning-controllability. 



IV. NUMERICAL RESULTS 

In this section, by following the approach presented 
above, we present numerical evidence of the useful- 
ness of R N+1 as an index for evaluating the pinning- 
controllability of a given complex network. Moreover, we 
show the behavior of R N+1 under variations of the con- 
trol gain re and the pinning probability p, in the case of 
scale-free complex networks. 
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FIG. 2: (Color online) Curve level sets of the natural logarithm of E, asymptotic value of the control error, as a function of 
both k and <r, for a BA network of 10 3 nodes, with average degree (k) = 2, p = 0.1. The left (right) panel shows the case of 
Combinatorial (Normalized) Laplacian. 



To validate our theoretical findings, we consider the 
case of a Barabasi- Albert (BA) scale- free network [ltij . 
characterized by a power-law degree distribution P(k) ~ 
k~ 3 (we have checked our results to be reproduced in 
a similar qualitative way for different kinds of networks 
and lattices, including real- world samples of networks), 
with N identical Rossler oscillators placed at the net- 
work vertices. Namely, the dynamics at each node 
i is described by the following vector field: f(xi) — 
f(xii,Xi2,x ia ) = (-x i2 - x i3 , xn + 0.165x i2 , 0.2 + 
{xn — 10)2^3). The output function h has been chosen, 
as in [19], to be h(x) — Hx, where H is the matrix: 
((1 0),(0 0),(0 1) ), indicating that the os- 
cillators are coupled through the variables xn and 2^3 

(; = i,2,...,at + i). 

Two cases have been considered: (1) the Combinato- 
rial Laplacian defined as Cij = Cji — —1 if i and j are 
connected (i ^ j), otherwise, and Cu — Hi Mi (corre- 
sponding to a symmetric network); and (2) the Normal- 
ized Laplacian defined as Cij = — l/fc, if i is connected 
to j (i ^ j), otherwise, and Cu = 1 Vi, which corre- 
sponds to an asymmetric network configuration (for more 
details, see [Hj]). 

For simplicity, we evaluate here the effect of n pinned 
nodes with gains K\ = K2 = ... = n n = K , on the struc- 
tural parameters R N+1 and M N+1 . For both the Combi- 
natorial and the Normalized Laplacians considered here, 
the spectra are real and thus M N+1 = 25]. The main 
results are shown in Fig. [1] for symmetric and asymmet- 
ric topologies. In both cases, the eigenratio R N+1 (in 
FigsQJa) and[T][c)) is characterized by a minimum around 
a specific value of the control gain k. This has immedi- 
ate relevance as it suggests that by appropriately tuning 
the control gain, it is possible to enhance the pinning- 
controllability of the network. Specifically, we observe 
that either a too large or too low value of k can reduce 



the network pinning-controllability. 

As shown in Figs UK a) and[ljc), when p is low (p = 0.1 
in the figure), R N+1 is higher than the eigenratio R 
associated with the corresponding uncontrolled network 
(represented by the continuous horizontal lines in FigfT]). 
Intuitively this is not surprising since, in order to control 
the network, we are requiring one more eigenvalue (than 
in the case of generic network synchronization [10|) to 
fall into the stability region of the MSF. 

Note that, for every choice of the coupling gain a, the 
eigenratio R N+1 may be conveniently varied by choosing 
an appropriate control gain n. At the same time, for 
every k, a bounded region of stability can be defined on 
the complex plane, and a timely choice of a is necessary 
in order to place all the a\Xi inside it. Thus, the stability 
of the reference evolution is sensible to both the values 
of the control gain and the coupling gain and values of 
k and <7, which are either too large or too small, may 
prevent the stability over the reference evolution of the 
network trajectories. 

It is worth noting here that there is a prominent differ- 
ence between the effects that varying a and k can have 
on the pinning-controllability of a given network. Specif- 
ically the width of S is determined essentially by the par- 
ticular shape of the MSF (which can be numerically com- 
puted, once the forms of the functions / and h are given) . 
On the other hand, the width of /C depends on the spec- 
tral properties of the matrix associated to the extended 
network, defined as proposed in this paper. Specifically, 
when the MSF is non-monotone, as in the case of the xz- 
coupled Rossler discussed in this paper, the width of the 
interval of the values of k, for which the stability of the 
reference evolution can be guaranteed, depends simply 
on the eigenratio R + , as shown in Fig. HJa)-(c). 

The asymptotic value of the control error E = 

JATJN Si=l It +AT W X iW) ^ s (t)\\dt, With ||x|| = + 
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FIG. 3: Left pictures: Combinatorial Laplacian (symmetric 
topology), (a) Behavior of the eigenratio /-Ln+i/^2 of the 
Laplacian spectrum as function of the control gain k for a 
BA network, with average degree (k) =2, a = 0.30, K = 10. 

(b) Control error at regime E as function of the control gain 
k under the same conditions as in the upper plot. 

Right pictures: Normalized Laplacian (asymmetric topology). 

(c) R N+1 vs. k for a BA network, with average degree (k) = 2, 
a — 2.8, K = 1.5. (d) E vs. n under the same conditions as 
in the upper plot. 

In all the plots the legend is as follows: squares (triangles) are 
used in the case of random (selective) pinning. In (a) and (c) 
the horizontal continuous lines represent the eigenratio R of 
the corresponding uncontrolled networks. 



1^2 1 + 1^3 1 j nas been computed under variations of the 
control gain k in Figs. D^b) and DJd) (respectively for 
fixed a = 0.30 and 2.8) and both k and a in Fig. [2] Note 
that interestingly the control error E in Fig. [5] behaves 
locally as a convex function of [k, <j\. 

Finally, the effect of a variable number of pinned nodes 
(randomly selected within the network) is shown in Fig. 
[3J Observe that, by increasing the number of reference 
nodes, differently from the case of variable control gains, 
we recover a monotone behavior in the eigenratio R + , 
as shown in Fig[3j^a). Interestingly, by increasing p, it 
is even possible to make R N+1 lower than R N , i.e. con- 
trolling the network becomes easier than synchronizing 
it. Also in Fig. [3JJb), the asymptotic value of the control 
error E is observed to decrease for increasing values of p. 

In Fig. [3J a comparison between random pinning 
(squares) and selective pinning (triangles) is also shown. 
In the case of selective pinning, the nodes have been cho- 
sen in the order of decreasing degree. We observe that 
selective pinning yields a lower eigenratio R N+1 as well 
as a lower control error E over different values of p, as 



shown separately in Figs. [3ja)-(b). This is in agreement 
with the previously reported results where a different ap- 
proach was used to evaluate the effectiveness of pinning 
control schemes @. Moreover, this shows the benefits of 
taking into account some network topological features in 
choosing some suitable nodes to pin. 

From a control design viewpoint, we observe that, 
while increasing the control gains may lead to a loss of 
pinning-controllability of the network, applying a larger 
number of controllers is always an effective strategy. On 
the other hand, the requirement of pinning a large num- 
ber of nodes can make the technique unfeasible in those 
networks where altering the dynamics of too many nodes 
can lead to a loss of functionality of the network itself, 
or is too costly and unpractical. 

The case of Normalized Laplacian is shown in Figs. 
|3jc)-(d), where the selective pinning is confirmed to per- 
form better than the random one. This indicates that 
high degree vertices continue to be better suited to con- 
trol the network than the others, even when the total 
strength of the interaction at each vertex is rescaled by 
its degree. 

Following these preliminary results, we believe other 
network properties such as degree correlation [l7j . clus- 
tering [2l[, centrality [2(| and community structure [22j 
can influence the effectiveness of pinning-control schemes. 
This represents the subject of future research activities. 



V. CONCLUSIONS 

In this paper we have presented a novel theoretical ap- 
proach to describing the controllability of networks un- 
der pinning control schemes. We have defined the novel 
concept of pinning-controllability and characterized the 
pinning-controllability of a given network in terms of the 
coupling gain, the control gain and the number of pinned 
nodes. We found that this property can be analyzed 
by investigating the synchronizability of an appropriately 
extended network. For instance, Fig. 1 shows that there 
are values of the control gain which render the eigenra- 
tio or the control error minimal. We wish to emphasize 
that information such as this could be used to under- 
stand why some physical networks possess certain values 
of the coupling and control parameters. Moreover, the 
methodology presented here can be an effective tool for 
the design of pinning control schemes of many biological 
and technological networks. 
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